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Abstract. We present pulse phase resolved spectra of the hypergiant high mass X-ray binary GX 301-2. We observed the 
source in 2001 October with RXTE continuously for a total on-source time of almost 200ksec. We model the continuum with 
both, a heavily absorbed partial covering model and a reflection model. In either case we find that the well known cyclotron 
resonant scattering feature (CRSF) at ~35 keV is ~ although present at all pulse phases - strongly variable over the pulse: the 
line position varies by 25% from 30keV in the fall of the secondary pulse to 38keV in the fall of the main pulse where it is 
deepest. The line variability implies that we are seeing regions of magnetic field strength varying between 3.4 X 10'^ G and 
4.2 X 10'- G. 
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1. Introduction 

GX 301-2 (4U 1223-62) is a High Mass X-ray Binary sys- 
tem consisting of a neutron star accreting from the strong 
stellar wind of the early type B-emission line star Wray 977. 
Due to difficulties in classifying Wray 977, the sy stem param- 
eters a re rather uncertain. In a detailed analysis, ' Kaper et'aP 
lll995h classified Wray 977 as a Bl lae-H hypergiant at a dis- 
tance of 5.3 kpc; we will adopt these values throughout this 
paper. If this classification is correct, Wray 977 is one of the 
most luminous and massive stars in our galaxy, with a lumi- 
nosity of 1.3 X 10^ Lq and a mass in excess of 48 M© (best fit 
~75 Mq). With one of the highest known wind mass loss rates 
in the Galaxy (M ~ 10"^ M© yr"'), the very slow wind (vk, = 
400 km s"^) can easily feed the neutron star with enough mate- 
rial to explain the observed X-ray luminosity of lO-'^ergs"'. 
We note that Koh et al.. (.1997.) doubt these extreme param- 
eters for Wray 977 as they are only barely compatible with 
their ana lysis, and argue that the earlier distance of 1.8 kpc 
jParkes et al..l980) and a spectral type of B2 lae is correct, with 
a correspondingly reduced luminosity of the Be primary. Even 
with these more moderate parameters, however, the system is 
unique. Note that the following analysis does not depend on the 
detailed choice of system parameters. 
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The neutron star is in an ecce ntric (e = 0.462) orbit with 
a period of 41.498 ± 0.002 days jKoh et alJll9'97t ISato et ail 
1986). Due to the violent wind accretion the source is ex- 
tremely variable - X-ray luminosit y changes of a factor of fiv e 
within one hour being common ( Rothschild & Soo'n3ll987l) . 
The underlying orbital variability, however, always follows 
the same pattern (Fig. 0: Shortly before periastron passage, 
the neutron star intercepts the gas stream from Wray 977 
( Leahv 1991, 2002), resulting in an extended X-ray flare during 
whic h the luminosity in creases by a factor of ~25 (see Fig.^ 
and iPravdo et aljfl99 5*). During periastron passage, the neu- 
tron star passes through the dense in ner wind of Wray 977, 
~0. 1 R* above the stellar sur face rsee lPravdo & Ghoshll2()0ll 
and Fig. 12}, and the X-ray activity decreases until it reaches a 
minimum after the periastron passage. This behavior is proba- 
bly due to the neutron star being almost eclipsed by Wray 977 
llLeahvll2002h and an intrinsically lower luminosity compared 
to the flare. Following this minimum, the X-ray luminosity in- 
creases slowly again over the orbit with a possible second max- 
imum near apastron, where the neutron star intercepts the spiral 
shaped gas stream from Wray 977 a second time ( Leahy 20()2). 

Using the Rossi X-ray Timing Explorer (RXTE) we ob- 
served GX 301-2 starting in 2000 October (JD 2451829.8 to 
JD 2451868.3) covering the greater part of the pre-periastron 
flare and the complete periastron passage in one contiguous 
~200ksec long observation (Fig. |3}. We used the same pro- 
cedures as in our earlier analysis of RXTE data from Vela X- 
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Fig. 1. Folded light curve using all available data (starting in 
1996 until mid 2003) on GX 301-2 of the All Sky Monitor 
on board RXTE. The light curve has been folded with the or- 
bital period of 41.498d (Koh et al. 1997). The peri astron pas- 
sage h as been extrapolated based on the ephemeris of lKoh et alJ 
lll997l) . The flare shortly before the periastron passage is very 
evident. Note the extended low following the periastron pas- 
sage which is probably due to the optical companion almost 
eclipsing the neutron star. For clarity the folded light curve is 
shown twice. 
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Fig. 2. Sketch of the system GX 301-2/Wray977 based on 
the parameters of Kaperet al. ( 1995). The neutron star passes 
Wray977 at a distance of < 0.1/?* during the periastron pas- 
sage. The time of the observation is marked by dashes and 
stronger line thickness, covering a significant part of the orbit 
due to the high velocity of the neutron star during periastron 
passage. The grey inner circle represents the size of Wray 977 
when usin g the old value s of Parkes et al.. C1980,) which are also 



when u smg the old value s < 
used bv lKoh et alJ ( Il997l) . 
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Fig. 3. X-ray Ught curve of the periastron passage of GX 301-2 
as observed with RXTE. The count rate is given in counts/sec 
per Proportional Counter Unit (PCU). Note the pre-periastron 
flare, the low following the periastron passage, and the over- 
all variability of the source. The time of the periastron passag e 
has been extrapolated using the ephemeris of lKoh et alJ ( ll997l) : 
the horizontal bar represents the uncertainty due to the uncer- 
tainties of the ephemeris. Our observation covered a significant 
part of the pre-periastron flare and the periastron passage itself. 
During the flare, count rates are up to 20 times higher than in 
the low following the periastron passage. 

Table 1. Systematic errors applied to the PCA-data to ac- 
count for uncertainties in the PCA calibration. We derived 
these values by fitting a two power law model simultaneously 
to a spectrum of a publ ic RXTE observation of th e Crab. See 
IWilms et al] ( 1 19991) and lKrevkenbohm et alJ ( l2002i) for a more 
detailed discussion of this procedure. 
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1 jRrevkenbohm et al] l2002h . Due to the improved calibration 
of the Proportional Counter Array (PCA), the systematic un- 
certainties associated with the PC A differ with respect to our 
earlier work (compare Table ]T\ and Krevkenbohm et al. '2002i 
Table 1). All analysis wa s done using FTOOLS 5.2 and XSPEC 
??.2.0aH llArnaucjl994 . Phase resolved spectra were produced 
with a modified version of the FTOOL fasebin, whic h prop- 
erly accounts for the HEXTE dead time jKrevkenbohm et alj 
2002). 

The remainder of this paper is structured as follows. In 
Sect. 12] we discuss the pulse period and pulse profiles. Sect.|3] 
discusses the different possible spectral models and detailed fits 
to phase resolved spectra with special emphasis on the behavior 
of the cyclotron fine. We summarize our results in Sect.|4| 



Ingo Kreykenbohm et al.: The variable cyclotron line in GX 301-2 



3 



2. Pulse profiles 

With a pulse period of ~700sec fewa nk et alJ [197^, 
GX 301-2 is one of the slowest pulsars kno wn. Since the firs t 
measurement of the period with Ariel-5 iWhite et alJll976l) . 
it has undergone dramatic variations: the period was stable at 
~700 sec until 1984, then the pulsar was spun up between 1985 
and 1990 jKoh et alJll997tlBildsten et alJll997i) to a period of 
~675 sec. In 1993, the spin behavior reversed (I Pravdo & GhoshI 
1200 ih and the pulsar has since been spinning down almost con- 
tinously: we determined the pulse period of our observation 
(JD 245 1830.8) to be 684.2 sec (see below). These long lasting 
spin-up and spin-down trends are in favor of GX 301-2 being 
predominantly a disk accreting system. Similar to other HMXB 
systems like Vela X-1, however, the pulse to pulse variations 
are very strong and the overall short term evolution of the pulse 
period is best descr ibed by a random walk model as is typical 
for a wind accretor jde Kool & Anzeiil993l) . 

To obtain pulse profiles, the bin (PCA) or photon arrival 
(HEXTE) times were first barycentered and then corrected 
for the orbital motion of GX 301-2 using the ephemeris of 
[Koh et al. ( 1997). We used epoch folding (Leahy et al. 1983) 
to determine the pulse period during the periastron passage to 
be 684.2 +0.2 sec, thus confirming that GX 301-2 is continuing 
its spin down. 

We used this period to create energy resolved pulse pro- 
files. The long period, the pulse to pulse variations, and the 
strong variability of the source exacerbate the determination of 
a pulse profile. In our case, GX 301-2 is brighter by a factor 
of more than ten during the pre-periastron flare compared to 
the low during the actual periastron passage (Fig.|3|l. We there- 
fore created several sets of energy resolved pulse profiles us- 
ing PCA -data: For the first set we used data taken during the 
flare only (see Fig. |3, while we used data from the low and 
directly after the pre-periastron flare for the second and third 
set (Fig.O. In general, the average pulse profile of GX 301-2 
is very well defined and similar to that of Vela X-1: it con- 
sists of two strong pulses at higher energies evolving into a 
more complex shape with substructures below 12keV. As one 
of the peaks is stronger in all energy bands than the other, we 
designate this peak as the main pulse and the weaker peak as 
the secondary pulse. The secondary peak is strongly energy- 
dependent: it is very weak (compared to the main pulse) at en- 
ergies below lOkeV. At energies above ~20keV it is signifi- 
cantly stronger, in fact, almost as strong as the main pulse. At 
these higher energies, the profile is approximately sinusoidal; 
however, the main pulse remains stronger throughout the en- 
ergy band covered here. Also similar to Vela X-1 is that the 
long term pulse profile is very stable. There are, however, sig- 
nificant short term variations. 

The interpretation of this energy dependence of the pulse 
profile is not straightforward. Above ~20keV, the sinusoidal 
pulse shape can be attributed to contributions from the two 
magnetic poles of the neutron star. It has been argued that 
anisotropic propagation of th e X-rays in the magnetized 
plasma of the accretion column llNageil98ll) or a complicated 
configuration of the magnetic field ("Mvtroph anov & Tsvg ^ 
Il978.) could produce the complex pulse profile at lower en- 
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Fig. 4. Pulse profiles in four energy bands from data taken dur- 
ing the pre-periastron flare (see Fig. E}- For clarity the pulse 
profiles are shown twice. Note that error bars are shown, but 
are too small to be seen in print. 

ergies as well. Alternatively, in analogy to Vela X-1, which 
also shows a strong energy dependence of the pulse profile 
(iKrevkenbohm et al. 1999, 2002), increased photoelectric ab- 
sorption or the accretion column passing through the line of 
sight can also be invoked to explain the observed profiles 
(Nagase et al. 1983). Due to the large and Compton thick ab- 
sorbing column during the flare (A^h ~ 10^^ cm"^, see be- 
low) it seems more likely that this latter explanation holds for 
GX 301-2 as well. To investigate the causes of the pulse profile 
variations, a detailed spectroscopic analysis of the pulse phase 
resolved X-ray spectrum is required. 

3. Phase resolved spectra 

3.1. Introduction 

As mentioned above, in the standard picture of accreting X- 
ray pulsars the strong magnetic field of the neutron star cou- 
ples to the infalling material and channels it onto the two 
magn etic poles where two hot spots emerge ( Burnard et alJ 
Il99ll) . If the magnetic axis is offset from the rotational 
axis, the ro tation of the neutron star gives rise to pulsations 
jPavidson & Ostriker 1973). The X-ray spectrum observed 
from the accretion column and the hot spot has generally been 
modeled as an exponentially cutoff power law, which is modi- 
fied by strongly varying photoelectric absorption from the ac- 
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Fig. 5. Definition of the six phase bins used for phase resolved 
spectroscopy using pulse profiles of data taken after the pre- 
periastron flare (see Fig. 0. Different shadings show the six 
phase bins in use: rise and fall of the main pulse (MPR and 
MPF), rise and fall of the secondary pulse (SPR and SPF), and 
the two pulse minima (OFFl and OFF2). 



cretion stream or the stellar wind at low energies (^ White et alJ 
[l980. J^983). Furth ermore, an Fe Kg fl uorescence line is ob- 
served at ~6.4keV. ISwank et alJ ( Il976l) also observed an iron 
edge near 7 keV. 

The strong magnetic field {B ~ lO'^ G) at the neutron stars' 
magnetic poles leads to the formation of cyclotron resonant 
scattering features (CRSFs, "cyclotron lines"). These lines re- 
sult from resonant scattering of electrons in La ndau levels in 
the ~10'^G magnetic field of neutron star (see ICoburn et all 
l2002t lArava&Hardinglll999t iMeszaros & Nage]lll985[ and 
references therein). The energy of the feature is nominally 
given by 

£c = ll-6keVx ^— X — ^ (1) 

where Ec is the centroid energy of the CRSF, z is the gravi- 
tational redshift at the scattering site (if the scattering occurs 
at the surface of the neutron star, z ~ 0.3, depending on the 
equation of state), and B the magnetic field strength. The en- 
ergy of the fundamental CRSF thus gives a direct measure of 
the magnetic field strength in the line forming region. 

Since the physical conditions are expected to vary over 
the emission region, the X-ray spectrum is expected to 
change with the viewing angle and therefore with pulse 
phase. This change is especially important when analyzing 
CR SFs since their shape depends strongly on the viewing an- 
le ("Arava-Gochez & Harding"200d'; 'Arava & Hai-ding"l999'; 



Isenberg et al. 1998a). thus making a spectrum averaged over 
the whole pulse period difficult to interpret. We therefore ana- 
lyzed the data from GX 301-2 using six phase intervals as de- 
fined in Fig.|5] the rise and fall of the main pulse, the rise and 
fall of the secondary pulse, and the two pulse minima. These 
intervals were chosen to give a good coverage of the pulse with 
a similar signal to noise ratio. 



3.2. The spectral model 

The lack of theoretical models of high enough quality for com- 
parison with observational data forces the use of empirical 
spectral models for the charac terization of the observed con- 
tinuum dKrevkenbohm et alJl999t) . Here, we describe the con- 
tinuum produced in the accretion column of the neutron star, 
^NS(£), by a power law with photon index F which is cutoff by 
the Fermi-Dirac cutoff" (iTanakai 19861) . 



^cont(£') - ^PL - 



(2) 

exp((£'-£'ciit)/£'F)+ 1 

where ApL is a normalization constant, and where Ecut and Ep 
are called the cutoff" and folding energy of the Fermi-Dirac cut- 
off". This continuum shape is then modified by a CRSF mod- 
eled as an absorptio n line with a Gaussian optical depth profile 
dCoburn et alJ2002l Eqs. 6 and 7), 



tgabs(£') = t/c X exp 



1 l E-Ec Y 
"2 1 ^rc ) 



(3) 



where Eq is the energy, crc the width, and dc the depth 
of the CRSF. The shape of this absorption line is sim- 
pler than the pseudo-Lorentzian line shape (as used by the 
XSPEC model CYCLABS and in earlier works; see, e.g., 
iMakishima & Ohashi|[l990l) and gives equally good fits. The 
overall X-ray spectrum emitted by the neutron star is then given 
by 

Ins(E) = /cont(£) exp (-TGABS(£)) (4) 

where I^oni is the continuum before scattering as given by, e.g., 
Eq.|2] Due to the strong stellar wind of Wray 977, however, the 
observed X-ray spectrum of GX 301-2 is strongly modified by 
photoelectric absorption, especially during times when the neu- 
tron star is close to Wray 977, such as in our observation. As 
the material enshrouding the neutron star is most likely clumpy, 
we expect (part of) the neutron star's X-rays to be strongly ab- 
sorbed close to the neutron star (A^h,2 in Eq.|5|i. Furthermore, all 
X-rays emerging from this zone will also be subject to photo- 
electric absorption in the overall stellar wind of Wray 977 (A^h,i 
in Eq.|5jl. A model describing such a physical situation is the 
absorbed partial covering model, which has the form 

/obs(£) = e-'^^^'^"' (l + ce-'^^'^"-^) /ns(£) + If.(E) (5) 

where crbf is the bound free absorption cross section per 
Hydrogen atom ( Wilms et al. 2000). A^h,i and A^h,2 are the col- 
umn densities of the two components, and c is the covering 
fraction of the absorber responsible for A^h,2- The possibility of 
a fluorescent Fe Ka line is accounted for by a Gaussian emis- 
sion line, /pe. 

A second possible scenario takes into account that during 
periastron passage the neutron star is extremely close to the sur- 
face of Wray 977. A large fraction of the neutron star's X-rays 
are thus intercepted by Wray 977 where they are backscattered 
by Thomson scattering or absorbed. This situation is analo- 
gous to th e Compton reflection hump observed in active galac- 
tic nuclei ("Lightman & White" 198 8*). In principle, the X-rays 
backscattered from within the atmosphere of Wray 977 could 
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be responsible for a significant part of the total observed X-ray 
flux. A model describing this physical situation has the form 

/obs(£) = e-""^'^" (/ns(£) + R/,s (E; n/2K) + h,iE)) (6) 

where all symbols are as described above and where 
R/(£'; Q./2n) describes the spectrum resulting from Compton 
reflecting X-rays with spectral sh ape I{E) of f gas with covering 
factor D./2n (Magdziarz & Zdzi arskilll995l) . Again, this spec- 
trum is photoabsorbed in the strong stellar wind of Wray 977. 

We note that these two pictures are not exclusive: the ac- 
cretion stream might well be patchy close to the neutron star, 
and the X-ray spectrum emerging from the accretion stream 
would then be reprocessed in the atmosphere of Wray 977. 
Detailed spectral fitting shows, however, that applying either 
of the above models already results in x^g^ ~ 1 • Combining 
the partial covering model and the reflection model is there- 
fore not necessary. With the available spectral data alone, it is 
not possible to test which of the two continuum models is the 
more realistic one - the absorption column of the partial cov- 
ering model is so high that it mimics the features of the soft 
end of the reflection spectrum (note that the shape of the re- 
flection spectrum below lOkeV is essentially identical to the 
shape of the incident continuum, absorbed by a column equiv- 
alent to one Thomson optical depth). Since the major aim of 
this work is the study of the cyclotron line, which is not af- 
fected by the details of the modeling of the continuum below 
lOkeV, we will concentrate our discussion in the remainder of 
this paper on the behavior of the source as reflected in the fits 
of the absorbed partial covering model (Table 1^. The corre- 
sponding fit parameters from the reflection model are listed in 
TablelH 

To allow the comparison of the spectral shape used here 
with other parameterizations of the continuum, we also mod- 
eled the data with several other standard pulsar continua (with 
and without a CRS F), especially using the high energy cutoff 
toite et al."l983, see Table 13. Residuals for these continua 
for the rise of the secondary pulse are shown in Fig.|6l In this 
pulse region (phase region SPR, see FigO the CRSF is very 
prominent such that the data are especially suited to illustrate 
the difficulties of the standard spectral models. 

3.3. Spectral fits 

Since the source is extremely variable, we accumulated phase 
resolved spectra not only for the whole data set, but separately 
also for the pre-periastron flare and the low following the flare. 
With a total usable exposure time of 176 ks for PC A and 89 ks 
live time for HEXTE the phase resolved spectra from the whole 
data set provide excellent statistical quality (note that the expo- 
sure time in HEXTE is typically lower by ~40% due to the 
dead time of the instrument; during our observation, the or- 
bit of RXTE was very close to the South Atlantic Anomaly 
thus giving rise to an increased particle background resulting 
in a deadtime of about ~50% in the HEXTE). In the following 
paragraphs, we will use this complete dataset for our analysis. 
See Tables|3]and|3for a list of the spectral parameters. 

The total usable exposure time during the flare is ~25 ksec 
for PC A and ~ 12 ksec for HEXTE. Even though this is much 
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Fig. 6. Residuals of fits of conventional neutron star continuum 
models to phase bin SPR (see Fig.|5}- Unless otherwise noted, 
the residuals shown include photoelectric absorption, an addi- 
tive Fe Kg line, and the CRSF a: Negative Positive Exponential 
(NPEX,|Maara 1995: M ihara et alJll998h . with strong residu- 
als especially below lOkeV. Despite the inclusion of a CRSF at 
32keV, significant residuals remain above 40keV (x^^^ = 17.2 
before, = 10.7 after including the line), b: power law with 
a Fermi-Dirac cutoff (FDC, Tanaka 1986|, x^^d ~ ^-^ without 
a CRSF (not shown), 5.5 with CRSF at ~35key ) c: power 
law with high energy cutoff after IWhite et alJ ( Il983l HEC; 
X^g^ = 25.4), exhibiting a very prominent and spurious fea- 
ture. Because of this spectral artifact, this model should not 
be used for cyclotron lin e searches ( Kretschmar et al J 1 1 997t 
iKrevkenbohm et alJll999 ): for comparison with earlier obser- 
vations the best fit values are shown in Table|2] d: the Smoothed 
High Energy Cutoff (SHEC), the HEC model smoothed by in- 
cluding a Gaussian absorption line (GABS) at the cutoff energy 
(Coburn et al. 2002), describes the data very well after apply- 
ing a CRSF at 35 keV: xi^ = 2.6 before including a CRSF 
(not shown), and 1.7 after the inclusion; see also Table |2j e: a 
standard partial covering model (PC, essentially this is Eq. |5l 
setting A^H.i = 0) still shows significant deviations below 8 keV: 
^red = 11-0 without a CRSF (not shown), = 5.4 after in- 
cluding a 33 keV CRSF, which are remedied by using the f : 
Absorbed Partial Covering (APC) model of Eq.|5] x^^d =4.1 
without a CRSF (not shown), and;^^^^^ - 0.8 after the inclusion 
of a CRSF at 34keV. g: Reflection Model of Eq.|51 Residuals 
are comparable to those of the APC; therefore, we use the APC 
model along with this reflection model to analyze the data. 
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Table 2. Parameters for spectral modeling of data from phase 
bin SPR (for d efinition, see Fig.lS} with a power law and a high 
energy cutoff (White et al. 1983, HEC) or a smoothed version 
of this cutoff ( Coburn et al..2002. SHEC) with and without the 
inclusion of a CRSF. At lower energies all three models are 
modified by photoelectric absorption. 



Since the data from the flare are very similar to the data 
from the complete set we only discuss the spectral parameters 
for the complete data set (unless mentioned otherwise). The 
data from the extended low following the pre-periastron flare 
(see Fig. ^ and Fig. |3} have a low statistical quality and are 
therefore not discussed in detail. 
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Fig. 7. a: Data and folded model of the rise of the secondary 
pulse. The model is an absorbed partial covering model (for 
discussion, see text), b: Residuals for the model without a 
CRSF and c: with a CRSF at 34.2+,^^ keV. The inset shows 
the pulse profile of GX 301-2 in the energy range from 5 keV 
to 20 keV. The marked region is the phase bin under discussion 
- the rise of the secondary pulse (see also Fig.|5|l- 



shorter than that of the whole observation, the flare spectra 
dominate over the rest of the observation due to the much 
higher flux during the flare. Since the low luminosity post pe- 
riastron data have a much lower statistical quality, the uncer- 
tainties of all parameters from this time are fairly large and the 
CRSF could not be detected in any phase region with high sig- 
nificance. 



3.3.1. The X-ray spectrum below 10keV 

As is expected from the physical picture outlined in Sect. 13.21 
the spectral parameters identified with the stellar wind - A^h.i, 
A^H,2, and the parameters of the Fe Ka line - do not depend 
significantly on the pulse phase as they originate very far away 
from the neutron star compared to the accretion column where 
the X-rays and the pulse forms. 

The high column density of the stellar wind as measured 
during the whole observation, A^h.i ~ 2 x lO^-'cm"^ indicates 
how deep the neutron star is embedded in the wind. With a 

n column density of as high as 2.3 x 10^"* cm"^, the second com- 
ponent of the partial covering model is close to being Compton 

n thick. As explained above, this large value is responsible for 

- our inability to differentiate between the partial covering model 
and the Compton reflection model. 

n In the low following the pre-periastron flare, the absorb- 

- ing columns are measured to be even higher The column den- 
sity for the "weakly" absorbed component, A^h.i, varies be- 

- tween 3.5 x 10^^ and 4.5 x 10^^ cm"^ and is thus twice as 
: large as during the flare, and also the column density A^h,2 is 
I increased: it varies between 4 x 10^"* cm"^ to over 10^^ cm"^. 

^ However, the observed count rate is not only due to the strong 

J photoelectric absorption but also due to an intrinsically lower 

i source luminosity: while the unabsorbed flux during the flare 

J is about 1 X 10"** ergs"', the absorbed flux, however, is ~ 

J 4.9 X 10"^ erg s"' . During the low the unabsorbed flux is lower 
by a factor of three (~ 3.5 x 10"^ erg s"') and the absorbed flux 

J by a factor of more then six (~ 7 x 10"'*' ergs s"'). 

J The high column density of the wind is also responsible 
for the strong Fe Ka fluorescence line. The line energy of 
6.48!QQ2keV, corresponding to an ionization level from Fev 
to Fe IX, implies that the iron, and thus probably all circumstel- 
lar material, is mildly ionized. The width of the line, crpe - 
0.33!Q Q2keV, is close to the energy resolution of the PCA, the 
very small uncertainties indicate that the line width is indeed 
measured. 



3.3.2. The pulsar continuum and the cyclotron line 

Pulse phase spectroscopy allows us to study the variation of 
the pulsar emission over the X-ray pulse. Our data analysis 
shows that the photon index F is almost constant over the pulse 
and only varies between -0.2 + 0.1 and 0.1 +0.1. Note that 
the spectrum is slightly but systematically softer in all phase 
bins during the pre-periastron flare. The parameters of the cut- 
off also do not vary over the pulse: The folding energy is 
only variable between 5.0 keV and 5.9 keV, while the cutoff en- 
ergy Ecm is consistent with a constant value of 15.5 keV ex- 
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cept for the rise of the main pulse where it is marginally lower 
(£cut = 12.9+!:^ keV). 

This continuum model alone, however, does not result in 
an acceptable fit in any of the phase bins (see e.g.. Fig. 03). In 
all phase bins a large absorption line like structure is present 
between 30keV and 40keV (see Fig.|5Jl. After the inclusion of 
a CRSF, the fits are acceptable in all phase bins with;^'^^^ ~ 1.0 
(see Tables |3] and 1^. The imp rovement is extremely signifi - 
cant: according to the F-Test ( Bevington & Robinsonl[l992i) ' 
the probability of a chance improvement is 2.5 x 10"" in the 
fall of the secondary pulse (the bin with the smallest improve- 
ment). We remark that the CRSF is found at almost the same 
energy and the same width (within uncertainties), but with a 
higher J when using the SHEC model applied in earlier pub- 
lications (compare Tables|5]and|3. 

Contrary to the continuum, the depth and the position of 
the CRSF (see Eq. |3 for definition of the CRSF parameters) 
are strongly variable with pulse phase (see Tables |3| E and 
Fig llOl l. This behavior is fairly typical of most CRSF sources 
(e.g.. Vela X-1 or Her X-1). We find that the line is deepest in 
the minimum between the main pulse and the secondary pulse 
(phase bin OFFl) with dc = 0.32^g }^, and it is also of com- 
parable depth in the adjacent phase bins, the fall of the main 
pulse (MPF, dc - 0.29;^^ ^^) and the rise of the secondary pulse 
(SPR, dc = 0.23+° JJ;^). While the depth of the CRSF in these 
three phase bins has a lower limit of 0.20, the CRSF is signifi- 
cantly less deep in the other three phase bins. At the same time, 
the position of the CRSF changes from Ec = 30.1+g7keV in 
the fall of the secondary pulse and Ec = 31 .0+S , keV at the 
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rise of the main pulse to significantly higher energies in the fall 
of the main pulse and the pulse minimum, where it becomes 
as high as ~ 37.9 keV (Fig. \IQI . Assuming a canonical mass 
of I.4M0 and a radius of 10km for the neutron star, the typ- 
ical gravitational redshift amounts to ~29% at the surface of 
the neutron star, such that the measured values imply a mag- 
netic field between B = 3.4 x 10'^ G and B = 4.2 x 10'^ G for 
GX 301-2. 

This change of the CRSF energy is correlated with its 
depth: the line is deepest in the phase bins where the CRSF 
is at the highest energies (correlation coefficient > 0.97). Note 
that this behavior is not due to the intrinsically lower signal 
during the interpulse: during the fall and the rise of the main 
pulse, where the data are of comparable statistical quality (and 
the source has a similar flux), we find a significant difference 
of the spectral parameters. We finally note that the line width 
of the CRSF crc is always in the range from 3 keV to 7 keV, 
typical for CRSFs at these energies. 

As we illustrate in Fig. |5J in addition to the changes of 
the depth and position of the CRSF, there are indications for 
a change in the shape of the feature itself. After fitting the 
CRSF in the rise of the main pulse some residuals remain (see 
Fig. Ell, indicating that the CRSF has a non-Gaussian line shape. 
Furthermore, the remaining residuals at ~50keV could indi- 



' Note that in even if systematic uncertainties are not an issue, how- 
ever, the extremely low false alarm probabilities make the detection of 
the line stable against even crude mistakes in the computation of the 
significance. 



Main Pulse Rise 




Lc) 



10 100 
Channel Energy [keV] 

Fig. 8. a Data and folded model of the rise of the main pulse for 
the absorbed partial covering model (see text for discussion), 
b Residuals for a fit without a CRSF, and c with a CRSF at 
3 1 .0!°^ keV. The inset shows the pulse profile of GX 301-2 in 
the energy range from 5keV to 20keV. The marked region is 
the rise of the main pulse. Note that fitting a Gaussian shaped 
CRSF does not completely remove the residuals of the CRSF; 
the resulting fit is acceptable, but not very good (x^^a - 1-2). 
See text for a discussion. 



cate the presence of an emission wing. The existence of such 
emission wings and overall variations of the line profile over 
the pulse have been pred icted by the numerical simulations of 
lArava & Harding! (il999l) . 

We also searched for any yet undetected second CRSF. As 
the addition of a second CRSF at any energy does not improve 
the fit, we tried to force an additional CRSF on the best fit at 
twice or half the energy of the ~35keV CRSF, especially in 
the phase bins where the x^^,^ is > 1-0 like the rise of the main 
pulse. A line between 15 keV and 20keV can be excluded with 
very high confidence: the upper limit for the depth of a line 
in this energy range is less than 10 A CRSF at twice the 
energy (between 60keV and 80keV) can not be excluded with 
such a high confidence, but it is still very unlikely that there 
is a secondary line present; the upper limit for the depth of 
such a line is 0.3, which given the low statistical quality of the 
data at 70keV is very low. Typical line depths in this regime 
are expected to be -1.0: a CRSF with a depth of 0.3 and a 
width of ~4keV can be fitted almost everywhere above 70keV 
without affecting the resulting significantly. Therefore, we 
conclude that no other CRSF is present in the data; we caution, 
however, that the source is only weakly detected above ~60keV 
and therefore we cannot exclude the presence of a feature in 
this energy range. 



Table 3. Fitted parameters from pulse phase spectroscopy for using the whole 200 ksec long observation of the periastron passage and the pre-periastron flare (see Fig.|3. The 
model is the absorbed partial covering model described in Sect. l3.2l bv Eq.|5]with and without the inclusion of CRSF. Note that the inclusion of the CRSF between 30keV and 
40keV in this data set improves the fits in all phase bins including in the pulse minima, where CRSFs are often either not present or insignificant (see e.g., iKrevkenbohm et all 
l2002h . All uncertainties quoted in this table and elsewhere in this paper are 90% confidence, dof are the degrees of freedom. The flux is absorbed flux given in erg s ' in the 
2-lOkeVband. 



Phase 


A'h,i[1022] 


A'h,2[1022] 


F 


Ecut [keV] 


£f [keV] 


Ef, [keV] 


o-Fe [keV] 


dc 


Ec [keV] 


<rc 




2 

^red 


Flux 


MPR 


22.7+2-3 
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-0 12+""'' 


q q+0.7 
^■^-1.3 
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3^+0.03 
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20.5!!;« 


218.8!!« 
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Table 4. As Tabled but for the reflection model. 



Phase 


A'h [10-2] 


F 


EcM [keV] 


Ep [keV] 


Refl 


£Fe [keV] 


ff-Fe [keV] 




Ec [keV] 
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Fig. 9. Variation of the CRSF for the six phase bins defined in 
Fig. Note that not only the depth and the energy of the CRSF 
changes with phase but also the actual shape of the line itself 
seems also to be variable. See text for a discussion of these 
issues. 
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Fig. 10. Variation of the energy and the depth of the CRSF over 
the pulse for the APC model. Although the values are slightly 
different for the REFL model, the variation of the parameters is 
very similar. For clarity the pulse is shown twice, a shows the 
PCA -countrate. Note that error bars are shown, but they are too 
small to be seen in print, b shows the variation of the energy of 
the CRSF over the pulse. Note that the energy variation over 
the pulse is definitely non-sinusoidal and therefore not due to 
simple angle dependence, c shows the variation of the depth of 
the CRSF over the pulse. See text for discussion. 



4. Summary and discussion 

4.1. CRSF Variability 

The major result of this paper is the observation of a varia- 
tion of the energy of the cyclotron line by almost 8 keV (cor- 
responding to about 25%) from 30.1^^^ keV in the fall of the 
secondary pulse to 37.9^J gkeV in fall of the main pulse. This 
large variation is sim ilar to that found in several ot her sources, 
such as 4U 1626-67 jHeindl & Chakrabartvll99qt) with a vari- 
ation of 2 5%, 4U 0115 -H63 with a valuation of 20% in the 
harmonic ( iHeindl et arAl999: Santangelo et al. 1999 ). Cen X-3 
with a variation of 30% jBurderi et alJl200(]|) . or Her X-1 with 
-25% (f Gruberetal]l200lL 119801) . On the other hand, the en- 
ergies of the CRSFs in Vela X-1 varie s only by ~10% with 
pulse phase ( iKrevkenbohm et al]l2002l) . Consistent with ear- 
lier interpretations, the strong energy variation of the cyclotron 
line in GX 301-2 argues that during different phases of the X- 
ray pulse, regions with different magnetic fields are observed. 
There are several plausible scenarios for such a suggestion. 

First consider a simple scenario that we observe an accre- 
tion column in a pure dipole field. In such a case the observed 
25% change of the cyclotron line energy would translate into 
a change in the height of the scattering region of ~5 km (or 



~50% of the neutron star radius). Such a change is rather un- 
likely, given that accretion colu mns are generally predicted to 
have heights of less than 1km jBurnard et al]ll99lt lieckeJ 
'1998'; 'Bro wn & Bildstenlll998 '). Apart from these theoretical 
considerations, a very tall accretion column would also imply a 
rather uncommon magnetic field configuration for GX 301-2: 
For a ~5 km high accretion column, the emission pattern is very 
likely a pure fan beam. To entirely block the lower parts of the 
accretion column such that the X-ray spectrum is dominated 
by the "top" of the column requires that the angle between the 
magnetic field axis and the spin axis, /3, is at least 40°, and 
larger if relativistic light bending is included. Detailed mod- 
els for the emission pattern of XRBs, which include relativistic 
light propagation (Kraus et al. 1996; Blum & Kraus 2000), as 
well as the statistical anal ysis of accreting neutron star pulse 
profiles teulik et alJl2003 ^. however, imply /3 <20°. We con- 
clude, therefore, that the accretion column of GX 301-2 is 
much lower than what would be expected if the variation of 
Ec were due to a pure dipole field geometry. 

A second possibility for the cyclotron line variability, 
which has been discussed, e.g., for Cen X-3, would be the 
presence of a very wide polar cap. For a dipole configura- 
tion, .Burderi et aL t,200Q.) show that an opening angle of the 
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order of 55° is required to produce a 30% variation of the line 
energy. Such a large polar cap, however, is unlikely to pro- 
duce the observed X-ray pulse shape - the pulse profile is ex- 
pected to be much broader in this case. Furthermore, as already 
noted by Burderi et al. ( 2000), such large polar caps would im- 
ply that the line energy varies symmetrically about the pulse 
maximum, which is not the case for any of the pulsars with a 
strong variation of the line energy. Finally, we also note that the 
strong magnetic field confines the accreting plasma, resulting 
in a much srnaller opening a ngle of the neutro n star's hot spot 
jLitwin et alJl200lUBeckerl l998: Basko & Sun vaevll973) . As 
a result we conclude that the extended polar cap scenario is an 
unlikely explanation for the large variation of the energy of the 
CRSF 

A third and final scenario for the pulse variation is that the 
observed X-rays come from an accretion mound, as suggested, 
e.g., by Burnard et al. ( 1991). In this case the observed varia- 
tion of the line energy is mainly due to the higher multipole 
components of the magnetic field in the mound and the X-ray 
emissivity profile is a combination of a fan and a pencil beam. 

Parenthetically, we note that slight variations of the line 
energy are also predicted for the case of homogeneous 
magnetic fields, as the location, depth, and shape of the 
CRSF are predicted to change with viewing an^le, i.e., 
pulse phase ( Meszaros & Nagel 1985; Isenberg et al. 'l998bi 
lArava & Harding 1999: Arava-Gochez & Ha rding 2000) . The 
strongest of these effects are the emission w ings of the funda- 
mental C RSF caused by "photon spawning" jArava & Harding! 
[l999: Arava-Gochez & Hardingn .2000) . where electrons ex- 
cited into a higher Landau level by resonant scattering of a 
photon, decay into the ground state by emitting photons at an 
energy which is roughly that of the fundamental CRSF. As a 
result the fundamental line is predicted to be shallower and 
have a more complex shape than the higher harmonics. When 
folded through the HEXTE response matrix, these variable 
complex line profiles can result in slight changes of the mea- 
sured line energy. These changes, however, are much smaller 
than the 25% energy variation seen here, and can thus be ex- 
cluded as the cause for the observed variation. Nevertheless, 
the inspection of the residuals of the pulse phase resolved fits 
indicates a possible change of the shape of the CRSF over the 
pulse (Fig.l^. Similar t o, e.g., 4U 01 15-H63 llHeindl etalJl999t 
ISantangelo et al.'l999'). the presence of an unresolved complex 
structure of the fundamental could be postulated to explain 
the not fully satisfactory values in some phase bins (e.g., 
X^^^~l-2 for bin MPR (see Fig.|8}, as compared to;^'^^^ < 1.0 
in the other phase bins). With the increasing energy resolution 
of modern gamma-ray detectors, it is foreseeable that such pro- 
files will be resolvable with newer instruments such as, e.g., the 
SPI instrument on INTEGRAL. 



4.2. Implications of the Line variability 

In recent years, several groups have searched for correlations 
between the parameters of the CRSF and the X-ray contin- 
uum in the hope of deducing information about the line and 
continuum formation process tMakishima et al...l999: .Coburni 



1200 IL and references therein). A major disadvantage of these 
searches, however, has been that they typically used pulse aver- 
aged data instead of pulse resolved spectroscopy, complicating 
the interpretation of the correlations found. 

One of the most interesting new correlations found in the 
large sample of X-ray pulsar s observ ed with R XTE and ana- 
lyzed in a uniform manner by'Coburn'l'ioOl) and'Cob urn eTaD 
(2002) has been between the relative width of the CRSF, crc/E, 
and its depth, dc- As noted by Coburn et al.l (120021) . part of 
the crc/£'cyc-'^c-plane is inaccessible to RXTE due to observa- 
tional constraints; however, our data all fall into the range that 
is well observable. As shown in Fig. \^ our phase resolved 
results for GX 301-2 confirm the overall correlation: deeper 
CRSFs are generally broader We note that this agreement be- 
tween the phase averaged and t he phase resol v ed dat a suggests 
that the correlations found by ICoburn et alJ (l2002l) from the 
phase averaged data are indeed real and not due to effects of 
the averaging. 

Having many data points for one object helps in under- 
standing the cTc/Ecyc-dc correlation in terms of a physical 
model. As we have outlined in the previous section, it is likely 
that the X-rays observed from GX 301-2 are produced in an 
accretion mound of moderate height at the magnetic poles of 
the neutron star Consider first the case that only one pole is 
visible. Over the X-ray pulse we view this pole under differ- 
ent viewing angles, 0, where is the angle between the line 
of sight and the magnetic field in the accretion column. For the 
high temperatures of the accretion column, the line width is ex- 
pe cted to be dominated by the rmal broadening. As shown, e.g., 
bv lMeszaros & Nagell ([1985), the anisotropic velocity field of 
the electrons in the accretion column leads to fractional line 
widths of 



erg 
Ec 



cos^ 



(7) 



where kT^ is the temperature of the electrons along the mag- 
netic field lines. Basic Comptonization theory suggests that kT^ 
can be estimat ed from the folding energy of the pulsar contin- 
uum, Ep (e.g., Burderi et al. 2000, and references therein). If 
we furthermore assume that the seed photons for the Compton 
upscattering in the accretion column are created throughout the 
volume of the accretion column, then detailed Monte Carlo 
simulations, e.g., of I senberg et alJ i 1998b) show that the depth 
dc of the CRSF is expected to be largest for a slab geometry 
when the line of sight is almost perpendicular to the direction 
of the magnetic field. For a uniform temperature accretion col- 
umn, where kT^. is constant, these slab models thus predict an 
anti-correlation between dc and crc/Ec, in contradiction to the 
observation. On the other hand, assuming a cylindrical geom- 
etry, the depth dc of the CRSF is expected to b e largest when 
the line of sight is parallel to the magnetic field ('isenberg et alJ 
[l998bj Fig. 19). This geometry therefore predicts a correlation 
between dc and crc/Ec, in good agreement with the observa- 
tion. 

We note, however, that even phase resolved data are diffi- 
cult to interpret. In the above discussion, we assumed that the 
observed data are dominated by emission from one homoge- 
neous emission region. In reality, this is not the case and the 
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Fig. 11. Fractional CRSF width crc/Ec versus the depth of the 
CRSF for several accreting neutron stars from RXTE data. 
Diamonds: values derived bv'Cobur n et alJ ll2002h from phase 
averaged spectra. Filled circles: values derived from phase re- 
solved spectra for GX 301-2 (this work). 



observed data are generally a mixture of contributions from 
both magnetic poles, which could influence the observed cor- 
relation. While it is difficult to disen tangle the contributions 
from the two poles dKraus et al.lll996i: iBlum & Krausl l2000i 
it seems likely that the physical conditions at both poles are 
slightly diff'erent. We would then expect the parameters of the 
X-ray continuum emitted by each pole to be dififerent, which 
would be reflected by changes in the observed continuum pa- 
rameters. In GX 301-2, however, the continuum parameters as 
exemplified by Ep are remarkably stable, suggesting either that 
only one magnetic pole contributes to the observed data, or that 
the poles have similar temperature. For this reason, it seems un- 
likely that a mixture of flux from the two poles causes the ob- 
served (Tc/dc correlation. Given the contradiction of predicted 
spectra of a slab geometry with the observations, the slab ge- 
ometry can either be ruled out or additional physical processes 
are operating. 

Finally, we note a second correlation found in our phase 
resolved data that i s also present in the phase average data 
llCoburn et al. 1'2002'): a coiTelation between the line width crc 
and the energy of the line Eq shown in Fig.[^ In terms of the 
simple cyclotron line broadening theory of Eq. Qthis correla- 
tion is rather unexpected since such a strong correlation (0.93 
for GX 301-2) is only possible if cos 6 does not vary apprecia- 
bly or 6 is close to 90°, where a relativistic treatment of Eqn.0 
is required and crc becomes independent of cos 6. A variation 
of 6 is required, however, since strong X-ray pulses are ob- 
served. For theta close to 90°, however, we would not expect to 
see strong pulsations, contrary to what is observed. We there- 
fore think that smaller values of theta are more likely. As an 
example, for 6 - 60°, the correlation constrains the viewing 
angle to vary by +6° around its mean value, tight even in the 
case of a pure fan beam. 
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Fig. 12. CRSF width crc versus CRSF Energy Ec for the same 
sample of accreting neutron stars as in Fig. ^2 Note the very 
strong correlation between Ec and crc for GX301-2 (correla- 
tion coefficient 0.93; filled circles). 

4.3. Summary 

To summarize, our analysis of the RXTE data of GX 301-2 
shows 

1. The continuum of GX 301-2 is well described by an ab- 
sorbed and partially covered pulsar continuum model. An 
alternative explanation for the X-ray spectrum is a reflected 
and absorbed pulsar continuum. Current data do not allow 
us to distinguish between these alternatives. 

2. GX 301-2 shows a strong cyclotron line at an energy of 
~35 keV. The line energy, depth, and width are variable 
with the pulse phase and do not depend on the chosen con- 
tinuum model. 

3. The variability of the CRSF energy as well as the pulse 
profile suggest that the X-rays from GX 301-2 originate in 
an accretion mound, and that the change in energy is not 
due to a change in height in a dipole field. 

4. The correlation between the relative width of the line and its 
depth is suggestive of changes in the angle between the line 
of sight and the magnetic field at the neutron star pole, in 
agreement with the standard paradigm for pulsating X-ray 
sources, but contradiction with the numerical simulations 
for the depth of the line for certain geometries. 
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